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Several  major  inverse  methods  (Stommel- Schott  method,  Wunsch  method,  and  Bernoulli 
method)  have  been  successfully  developed  to  quantitatively  estimate  the  geostrophic 
velocity  at  the  reference  level  from  hydrographic  data.  No  matter  the  different  appear¬ 
ance,  they  are  based  on  the  same  dynamical  sophistication:  geostrophy,  hydrostatic,  and 
potential  density  (p)  conservation  (Davis,  1978).  The  current  inverse  methods  are  all 
based  on  two  conservation  principles:  potential  density  and  potential  vorticity  ( q  =  fdp/ 
dz)  and  require  /5-turning.  Thus,  two  necessary  conditions  can  be  incorporated  into  any 
inverse  methods:  (1)  non-coincidence  of  potential  density  and  potential  vorticity  surfaces 
and  (2)  existence  of  vertical  turning  of  the  velocity  (/5-turning.)  This  can  be  done  using  the 
P-Vector,  a  unit  vector  in  the  direction  of  VpxVq  (Chu,  1994, 1995).  The  first  necessary 
condition  becomes  the  existence  of  the  P-vector,  and  the  second  necessary  condition  leads 
to  the  existence  of  the  P-vector  turning  in  the  water  column.  Along  this  line,  we  developed 
the  P-vector  inverse  method  with  a  pre- requirement  check-up.  The  method  was  verified 
in  this  study  using  the  Modular  Ocean  Model  (MOM)  from  Pacanowski  et  al.  (1991) 
version  of  Bryan-Cox-Semtner  ocean  general  circulation  model  (OGCM),  which  is  based 
on  the  work  of  Bryan  (1969).  The  statistically  steady  solutions  of  temperature  and  salinity 
from  MOM  are  used  as  a  “no-error  data”  set  for  computing  absolute  geostrophic 
velocities  by  the  P-vector  inverse  method.  Circulations  are  similar  between  the  MOM 
statistically  steady  solutions  and  the  P-vector  solutions.  Furthermore,  the  quantitative 
analysis  shows  that  this  inverse  method  has  capability  of  picking  up  the  major  signal  of 
the  velocity  field. 
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1.  Introduction 

Our  understanding  of  the  mid-latitude  large-scale  ocean 
circulation  has  been  greatly  benefitted  by  a  remarkable  set 
of  papers  by  Stommel  and  collaborators  (Stommel  and 
Schott,  1977;  Schott  and  Stommel,  1978,  Beringer  and 
Stommel,  1980),  Wunsch  and  collaborators  (Wunsch,  1978; 
Wunsch  and  Grant  1982),  and  Killworth  (1986).  Their  work 
makes  it  possible  to  obtain  ocean  general  circulations  from 
observations  of  temperature  ( T)  and  salinity  ( S ).  The 
physical  base  for  calculating  geostrophic  velocity  from 
hydrographic  data  is  the  thermal  wind  relation 


u  =  u0 
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f  #*'■ 


fPo  z°  dy 
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V  =  v0 


fpQ  z°  dx 
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where  (u,  v),  (no,  vo)  are  the  geostrophic  velocity  at  any  depth 
z  and  at  a  reference  depth  zo,  p  is  the  in  situ  water  density, 
Po  is  the  characteristic  value  of  the  density,  and  /  is  the 
Coriolis  parameter,  which  is  a  function  of  latitude.  Here  the 
Boussinesq  approximation  has  been  used.  As  mentioned  by 
Olbers  etal.  (1985),  the  quantities  T,  S  are  relatively  easy  to 
measure,  and  in  contrast  to  velocity  observations,  the  cli¬ 
matological  signal  in  the  T,  S  fields  is  less  contaminated  by 
energetic  smaller-scale  motions  induced  by  eddies  and  waves. 
Equations  (1)  and  (2)  indicate  that  the  hydrographic  data 
only  determine  the  baroclinic  geostrophic  currents.  The 
reference  velocity  (wo,  vo)  still  needs  to  be  determined. 

Based  on  the  geostrophy,  hydrostatic  balance,  and 
mass  conservation,  several  major  inverse  techniques,  i.e., 
the  /5-spiral  method  (Stommel  and  Schott,  1977;  Schott  and 
Stommel,  1978),  the  Wunsch  method  (Wunsch,  1978),  and 
Bernoulli  method  (Killworth,  1986)  have  been  successfully 
developed  to  quantify  the  geostrophic  velocity  at  the  reference 
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level  (uo,  i;o)-  We  refer  reader  to  an  excellent  review  paper 
on  /5-spiral  method  by  Olbers  et  al.  (1985).  Davis  (1978) 
pointed  out  that  the  /5-spiral  method  and  the  Wunsch 
method,  no  matter  how  different  in  appearance,  are  based  on 
the  same  order  of  dynamical  sophistication  and  differ  from 
implicit  assumptions  about  the  scales  of  oceanic  variability 
and  different  definitions  of  the  smooth  field  to  which  the 
dynamical  model  pertains. 

As  pointed  out  by  Wunsch  and  Grant  (1982),  in  de¬ 
termining  large-scale  circulation  from  hydrographic  data, 
we  can  be  reasonably  confident  on  the  assumptions  of 
geostrophic  balance,  mass  conservation,  and  no  major  cross- 
isopycnal  mixing  (except  water  masses  are  in  contact  with 
the  atmosphere).  The  density  of  each  fluid  element  would  be 
conserved,  i.e., 


V  •  Vp  =  0 


(3) 


VpxVq *  0. 

Stommel  and  Schott  (1977)  pointed  out  that  because 
the  horizontal  component  of  velocity  rotates  with  depth  in 
the  open  ocean  (/5-spiral),  absolute  velocities  can  be  ob¬ 
tained  from  observations  of  the  density  field  alone.  Since  all 
the  inverse  methods  are  dynamically  equivalent,  the  /5-spiral 
should  be  satisfied  before  using  any  inverse  method.  This 
leads  to  the  second  necessary  condition. 

Necessary  Condition  2:  The  velocity  (u,  v)  should  have 
vertical  turning  (i.e.,  the  /5-spiral)  in  the  water  column,  i.e., 
the  horizontal  velocity  should  change  direction  with  depth 


u 


(k) 


u 


(m) 


(*) 

(m) 


(7) 


where  p  is  the  potential  density.  The  conservation  of  potential 
vorticity  equation  (Pedlosky,  1986)  can  be  obtained  by 
differentiating  (3)  with  respect  to  z,  using  the  geostrophic 
and  hydrostatic  balance,  and  including  the  latitudinal 
variation  of  the  Coriolis  parameter, 

V-V9  =  0.  q  =  f^-.  (4) 

oz 

Use  of f dpi dz  may  induce  a  small  but  systematic  error  into 
estimation  of  potential  vorticity  (Needier,  1986).  Equations 
(3)  and  (4)  indicate  that  V  is  perpendicular  to  both  Vp  and 
Vp,  therefore,  the  velocity  V  is  parallel  to  Vp  x  Vp. 

Recently,  Chu  (1994,  1995)  proposed  using  a  unit 
vector  (P-vector) 

p (5) 
Vp  X  Vq\ 


somewhere  in  the  water  column.  For  a  given  level  z  =  Zk,  if 
we  cannot  find  a  level  zm  such  that  (7)  is  satisfied,  the  inverse 
method  will  fail  to  get  velocity  at  Zk  of  this  water  column. 

Before  using  any  inverse  method,  we  need  to  check  if 
these  two  necessary  conditions  are  satisfied.  If  one  of  them 
is  not  satisfied,  we  cannot  use  any  inverse  method  to  obtain 
the  velocity  field  from  the  T,  S  fields  for  that  water  column. 
The  next  section  shows  that  a  recently  proposed  P-vector 
concept  (Chu,  1994,  1995)  can  bring  the  two  necessary 
conditions  into  the  /5-spiral  method. 

3.  P-Vector 

Existence  of  a  P-vector  (5)  implies  the  satisfaction  of 
the  Necessary  Condition  1.  This  provides  the  first  check 
point  to  see  if  the  potential  density  surface  coincides  with  the 
potential  vorticity  surface.  If  the  necessary  condition  1  is 
satisfied,  the  P-vector  lies  on  the  intersection  of  the  potential 
density  and  potential  vorticity  surfaces  (Fig.  1). 

The  relationship  between  velocity,  V  =  (u,  v,  w),  and 
P  =(PX,  Py,  PZ),  IS 


to  obtaining  North  Atlantic  ocean  circulation  from  hydro- 
graphic  data.  The  results  were  quite  agreeable  with  obser¬ 
vations.  The  major  purposes  of  this  paper  are  to  demonstrate 
the  benefit  of  using  P-vector  in  /5- spiral  method  and  to 
evaluate  the  P-vector  method  using  an  ocean  general  cir¬ 
culation  model. 

2.  Necessary  Conditions  for  the  /5-Spiral  Method 

The  three-dimensional  velocity  field  can  be  determined 
from  the  density  field  unless  the  potential  vorticity  and 
density  surfaces  coincide  (Stommel  and  Schott,  1977).  This 
leads  to  the  first  necessary  condition  for  the  validity  of  any 
inverse  methods. 

Necessary  Condition  1:  The  potential  density  surface 
does  not  coincide  with  the  potential  vorticity  surface;  i.e., 


V  =  r(x,y,z)  P 


(8) 


where  r  is  the  proportionality.  Applying  the  thermal  wind 
relation  to  any  two  different  depths  Zk  and  zm,  shown  as  in 
Fig.  2,  a  set  of  algebraic  equations  for  determining  the 
parameter  r  is  obtained 


M)p{k)  _  Jm)  p(m) 

1 X  ' 


An 


km 
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(k)  p(k)  (m)  p{m)  _  * 

1  ry  1  ry  ~  m  km 

which  are  two  linear  algebraic  equations  for  r(k>  and  rim).  Here 
r(l)  =  r(x,  y,  Zi ),  and 
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(12) 


{Aukm’Avkm)  =  -^-  JT 


d£_  _dp_ 

fPo  iz'n  {  dy  ’  dx 
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(10) 


As  soon  as  dk)  is  obtained,  the  velocity  field  can  be  com¬ 
puted  by  (8). 

If  the  determinant 
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the  algebraic  equations  (9)  have  definite  solutions  for  r{k) 
(m  ^  k ): 


Thus,  the  P-vector  has  the  second  good  feature  which 
is  the  check-up  for  the  /5-turning.  Consider  the  P-vector  at 
two  different  levels,  z  =  Zk  and  z  =  zm  (Fig.  2(a))  and  use  P/„ 
Pv  to  denote  the  horizontal  and  vertical  components  of  the  P- 
vector.  The  determinant  (11)  is  the  sine  of  the  vertical 
turning  angle  between  Ph(k)  and  Ph{m)  (Fig.  2(b)),  i.e., 


Fig.  1.  The  absolute  velocity  and  the  intersection  of  the  surfaces 
of  potential  density  and  potential  vorticity. 
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where  akm  indicates  the  /5-spiral  turning  angle  between  the 
two  levels  Zk  and  zm- 

In  order  to  use  any  inverse  method,  we  should  check  if 
the  inequality  (11)  holds.  If  (11)  is  not  satisfied,  i.e.,  the 
horizontal  velocity  does  not  have  /5-turning, 

sin  akm  =  0,  or  P^  =  const,  P^  =  const, 

for  all  m  (14) 

all  the  inverse  methods  fail.  This  is  the  second  necessary 
condition  for  checking  the  /5-turning. 


Y 


Fig.  2.  Vertical  turning  of  the  P-vector:  (a)  P  vector  at  two  different  levels,  and  (b)  turning  angle  between  two  levels  (from  Chu,  1995). 


Evaluation  of  P-Vector  Method  1 87 


4.  P-Vector  Inverse  Method 

For  water  columns  surviving  the  two  necessary  con¬ 
ditions,  we  may  use  ( 1 2)  to  compute  r(fe)  for  the  level  Zk-  There 
are  (N  -  1)  sets  (m  =  1,  2,  k  -  1,  k  +  1,  N)  of  Eq.  (9)  for 
calculating  r(k\  Here,  N  is  the  total  vertical  levels  of  the 
water  column.  All  the  (N  -  1)  sets  of  equations  are  com¬ 
patible  under  the  thermal  wind  constraint  and  should  pro¬ 
vide  the  same  solution.  However,  due  to  errors  in  mea¬ 
surements  (instrumentation  errors)  and  computations 
(truncation  errors),  the  parameter  r(k)  may  vary  with  m.  We 
have  developed  an  optimization  scheme  to  minimize  errors. 

If  the  absolute  velocity  (Ak\  v{k))  is  known,  we  may  use 
the  thermal  wind  relation  ( 1 0)  to  obtain  the  absolute  velocity 
at  any  level  m, 


Fig.  3.  Error  caused  by  uncertainty  of  the  velocity  at  the  level  k. 
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Table  1.  Distribution  of  vertical  levels. 


^m)  =  u{,n)dp(m)/dx  +  v{m)dp{m)/dy  (15) 
dp{m)/dz 

The  computed  [u^n\  v^m\  w^]  may  not  be  in  the  same 
direction  as  the  P-vector  [Px(mK  Pym\  PSm)]  at  the  level  m.  If 
we  assume  that  at  the  level  m,  the  P-vector  exists 


\VpxVq\ 

|vP|v?|  1 


(16) 


and  that  the  velocity  should  parallel  the  vector  P(m\  an  error 
can  be  easily  defined  by  (Fig.  3) 


Level 

Depth  (m)  of  T,  S 

Thickness  (m) 

1 

12.5 

25.0 

2 

37.5 

25.0 

3 

70.0 

40.0 

4 

125.0 

70.0 

5 

215.0 

110.0 

6 

370.0 

200.0 

7 

635.0 

330.0 

8 

1025.0 

450.0 

9 

1575.0 

650.0 

10 

2350.0 

900.0 

11 

3250.0 

900.0 

12 

4150.0 

900.0 

y(m)p(m)  _  y(m) 

=  2  sin 

~  i 

xT" 

3^ 

<  > 
3^ 

_ _ l 

L2 v  'J 

-V(w)sin(p(w),V(m)), 


Substitution  of  (17)  into  (18)  leads  to  a  2  x  2  algebraic 
(17)  equations  for  determining  ( u ('k\  v®), 


here  V(m)  =  Iv^l,  and  (P(m),  V^)  is  the  angle  between  two 
vectors  P(m)  and  The  total  error  of  the  water  column 
velocity  caused  by  the  uncertainty  of  ( u(k) ,  v®)  is 


Anu^  +  A12v^  =  Fx 
A2luw  +  A22vW  =  F2 


(19) 


E  = 


)2 


where  hm  is  the  thickness  of  the  m-th  layer  (see  Table  1). 

We  determine  the  velocity  (iEk\  v(fc))  such  that  the  total 
error  E  becomes  minimum,  i.e., 


c)E 

du(k) 


=  0, 


dE 

dv{k) 


0. 


(18) 


where 
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and 
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dp(m)  tdx 
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(21) 


The  absolute  velocity  at  the  level  z  =  Zk  can  be  computed 
by 


u(k)  _  FXA12  -  F2AX2  v(k)  _  F2An  -  FXA2X  ^2) 
F\  i 422  ~~  FnA2X  FnA22  -  Fl2A2i 

Substitution  of  (21)  into  (20)  makes  the  second  necessary 
condition  (14)  into  that  if 


5.  Modular  Ocean  Model  (MOM) 

Any  inverse  method  involves  two  different  kinds  of 
errors:  observational  and  modeling  errors.  The  best  way  to 
verify  the  model  is  to  use  a  no-error  data  set.  Since  there  is 
no  such  data  set,  we  may  use  a  set  of  steady  state  solutions 
from  a  numerical  model  as  a  no-error  data  set.  In  this  study, 
we  use  the  steady-state  solutions  of  temperature  and  salinity 
from  Pacanowski  et  al.  (1991)  version  of  the  Bryan-Cox- 
Semtner  ocean  general  circulation  model  (OGCM),  which  is 
based  on  the  work  of  Bryan  (1987).  The  model  domain 
consists  of  a  60°  square  box  in  latitude-longitude  space  from 
10°-70°N  and  10°-75°W.  Along  the  western  boundary,  an 
idealized  shelf  with  a  structure  similar  to  that  in  Holland 
(1973)  is  included.  The  horizontal  grid  spacing  is  2°  latitude 
by  2°  longitude.  The  model  has  12  levels  in  the  vertical,  and 
the  depth  distribution  is  the  same  as  that  of  Cai  (1995)  and 
are  listed  in  Table  1.  Values  assigned  to  the  various  model 
parameters  are  listed  in  Table  2.  No-slip  and  insulating 


AhA22  A12A21  -  0  (23) 


all  the  inverse  methods  fails,  In  other  words,  the  second 
necessary  condition  can  be  written  as 

|  A  1^22  ~  A12  Al  |  >  £2  (24) 

there  exists  /5-turning  and  we  may  use  inverse  method  to 
obtain  velocity  fields. 

The  P- vector  inverse  method  is  essentially  the  /5-spiral 
method  with  two  necessary  conditions.  The  benefit  of  using 
the  P-vector  method  is  to  filter  out  those  data  points. 


Table  2.  Model  coefficients. 


Parameter 

Symbol 

Value 

Horizontal  diffusivity 

Am 

1  x  10  3  m2s-1 

Horizontal  viscosity 

Amh 

2  x  10  5  m2s_1 

Vertical  diffusivity 

Ajy 

1  x  10"4  mV1 

Vertical  viscosity 

4  MV 

1  x  10  ^  m2s_1 

Characteristic  density 

Po 

1025  kg  m  ~3 

Heat  capacity 

r 

4  x  10  3  J  kg  ^K"1 

Fig.  4.  Surface  boundary  conditions.  Temperature  and  salinity  profiles  used  in  Haney-type  restoring  forcing  condition  and  surface  wind 
stress. 
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boundary  conditions  are  applied  at  the  lateral  boundaries. 
The  model  uses  the  Cox  (1987)  parameterization  to  compute 
vertical  diffusion  and  convection  implicitly.  The  enhanced 
vertical  diffusivity  in  regions  of  static  instability  is  set  at  105 
cm2s_1,  which  is  the  convective  adjustment  in  the  model. 

The  temperature  and  salinity  of  upper-most  level  are 
relax  under  a  Haney  (1971)  restoration  to  a  zonally  uniform 
temperature  profile: 

T  l(p)  =  30.0-^- (p  + 4.0 5.0)]  (25) 

v  '  80.0  L  75.0  v  ’] 


and  a  salinity  profile: 


Fig.  5.  The  x-dependent  bottom  topography  used  in  the  OGCM. 


Potential  Den&ly 


Fig.  6.  Statistically  steady  potential  density  (kg/m3)  from  the  OGCM  simulation  at  different  depths:  (a)  12.5  m,  (b)  37.5  m,  (c)  70  m, 
(d)  215  m,  (e)  1,575  m,  and  (f)  3,250  m. 


190  P.  C.  Chu  etal. 


Y  32  r  5  283  1  surface  years  (6,000  years  at  the  bottom)  the  model  reaches 

S{(p )  =  35.0-  (p  +  0.84 sin  ———((p- 15.0)  (26)  a  statistically  steady  state  (total  kinetic  energy  reaches 

w/  vy  •  «_/  *_/  •  W  •  -|  •  i  •  \ 

equilibrium). 

with  a  relaxation  time-scale  of  25  days  (over  a  modeled  top  MOM  Model  Generated  Steady-State  Data 

mixed  25  m).  The  model  is  also  subject  to  the  zonal  wind  The  M0M  model  out  t  includes  tential  te  a. 

stress  of  Bryan  (1987).  Figure  4  shows  the  surface  boundary  tuR,  m  salinity  (5)>  and  velocity  („_  w)  The  statistically 

conditions.  The  bottom  topography  is  assumed  only  depen-  steady  potential  denslty  field  was  computed  from  g  and  5 

dent  on  .r,  and  has  a  staircase-type  change  west  of  64°W  and  fields  Figure  6  shows  the  ^  (  _  1000)  flelds  at  different 

no  change  east  of  64  W  (Fig.  5).  depths.  The  most  evident  features  are  listed  as  follows.  The 

The  techniques  of  Bryan  ( 1987)  for  acceleration  equili-  maximum  values  (29.2-29.3  kg/m^)  appear  at  the  high 

bration  of  the  model  solution  are  used.  This  includes  using  latitudes(near70oN)atalldepths.Thehorizontal<j8gradient 

a  longer  time  step  at  depth.  The  acceleration  factor  increases  rediices  as  ,he  depths  increases  The  low  a$  water  appears  at 

from  1  at  the  surface  level  to  8  at  the  bottom  level.  After  750  the  equatorial  region  near  surface  and  stretches  northward  as 


Pdenlial  Vorticity  (1  .e-S  kg/(mfl  s)) 
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Fig.  7.  Statistically  steady  potential  vorticity  from  the  OGCM  simulation  at  different  depths:  (a)  12.5  m,  (b)  37.5  m,  (c)  70  m, 
(d)  215  m,  (e)  1,575  m,  and  (f)  3,250  m. 
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the  depth  increases.  Near  the  surface  (Fig.  6(a))  the  o q 
curves  are  almost  zonal  except  at  the  western  boundary, 
where  the  Oq  curves  bend  towards  the  north,  indicating  the 
northward  movement  of  low  Oq  equatorial  water.  At  the 
depths  of  37.5  m  (Fig.  6(b))  and  70  m  (Fig.  6(c)),  the 
northward  movement  of  the  equatorial  low  Oq  water  be¬ 
comes  more  evident.  At  the  depth  of  215  m  (Fig.  6(d)),  the 
low  Oq  zone  shifts  northward.  Two  centers  of  low  Oq  (27.4 
kg/m3)  are  found  at  the  southwest  corner  (25°N,  74°W)  and 
at  the  middle  of  the  eastern  boundary.  At  the  deep  levels 
(Figs.  6(e)  and  6(f)),  Oq  is  very  uniform.  The  minimum  zone 
of  Oq  is  located  at  around  60°N. 

The  statistically  steady  potential  vorticity  field  ( q )  was 
computed  from  Oq.  Figure  7  shows  the  q  ( =fdoQ/dz )  field  at 
different  depths.  The  most  evident  features  are  listed  as 
follows.  The  values  of  q  are  either  negative  or  zero.  There  is 


an  anticyclonic  gyre  occupying  the  majority  of  the  domain. 
In  upper  levels  (12.5  m,  37.5  m,  and  70  m),  the  strongest 
negative  q  center  (with  a  value  around  -1.4  x  10-6  kg  nr4 
s-1)  is  located  near  the  southwest  corner  (60°-70°W,  20°- 
35°N).  As  depth  increases,  the  magnitude  of  ^-values  re¬ 
duces.  The  strongest  negative  q  center  with  values  around 
-2.5  x  10~8  kg  m^s-1  at  1,575  m  depth  and  -2.2  x  10~9  kg 
m~4s_1  at  3,250  m  depth  (2-3  orders  of  magnitude  smaller 
than  the  upper  levels),  is  located  near  the  northeast  corner 
(20°-10°W,  55°-65°N). 

We  may  verify  the  MOM  p,  q  data  in  terms  of  the  two 
necessary  conditions.  The  vector  product  of  Vp  and  Wq  is 
written  by 

Vp  x  Vq  =P  IVpl  IVpIsinr) 


o 

— p - p - ■  ■  ■ - 

|  k 

r-. 

s 

- 

- 

■ 

■ 

„  o 

aLri 

- 

3 

■rd  rj- 

r 

- 

■ 

. 

8 

- 

■ 

■ 

o 

CM 

- 

-  4 

- 

■ 

- 

O 

- * - ■■  i _ i _ _ _ i _ 

oi — i —  ■  -  ■  ■  i  i _ . _  ■  -  i 

-70  -60  -50  -40  -30  -1C  -70  -60  -50  -40  -3Q  -20  -10 

Longitude  Longitude 


Fig.  8.  Regions  with  Isindl  <  10  5  at  (a)  370-m,  (b)  635-m,  (c)  1,025-m,  (d)  1,575-m,  (e)  2,350-m,  and  (f)  3,250-m. 
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where  8  is  the  intersection  angle.  Since  8  cannot  be  exactly 
zero  when  computed  from  data.  A  small  value  was  taken  as 
the  criterion  8C  =  10  5.  If  I  <51  <  8C,  the  p  surface  is  thought  to 
be  parallel  to  the  q  surface.  Figure  8  shows  the  regions  at  six 
different  depths  where  the  first  necessary  condition  fails. 

7.  MOM  Generated  Statistically  Steady-State  Veloc¬ 
ity  Field 

Figure  9  shows  the  MOM  generated  statistically  steady 
horizontal  velocities  at  several  depths:  (a)  12.5  m,  (b)  37.5 
m,  (c)  70  m,  (d)  215  m,  (e)  1575  m,  (f)  3250  m.  The 
circulation  patterns  can  be  outlined  as  follows.  Westward- 
moving  equatorial  currents  are  evident  in  upper  levels  (z  > 
-215  m)  and  disappear  in  the  deeper  layer  (z  <  -2 1 5  m) .  The 


width  of  the  upper  layer  equatorial  current  is  around  6°  in 
latitude.  When  the  upper  layer  equatorial  current  approaches 
the  western  boundary,  it  turns  direction  and  becomes  the 
western  boundary  current.  The  maximum  velocity  of  the 
western  boundary  current  reaches  1 3  cm/s.  The  intermediate 
level  (1,575  m)  and  deep  level  (3,250  m)  feature  a  westward 
current  originating  at  the  northeast  corner.  This  current  turns 
southwestward  at  around  50°W  and  branches  into  two 
western  boundary  currents  (northward  and  southward)  as 
approaching  the  western  boundary.  The  bifurcation  is  found 
at  40°N  for  the  depth  of  1,575  m,  and  at  60°N  for  the  depth 
of  3,250  m.  The  deep  currents  are  weak  (<2  cm/s). 

The  velocity  spiral  can  be  identified  for  each  grid  point 
by  the  plot  of  Vi,  of  all  the  depths  (Fig.  10).  There  are  12 
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Fig.  9.  Statistically  steady  velocity  from  the  OGCM  simulation  at  different  depths:  (a)  12.5  m,  (b)  37.5  m,  (c)  70  m,  (d)  215  m,  (e)  1,575 
m,  and  (f)  3,250  m. 
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Fig.  10.  Velocity  spirals  of  Vi,  (12  vectors  at  each  grid).  Two  numbers  given  at  each  grid  indicate  the  maximum  turning  between  two 
consecutive  levels. 


vectors  at  each  grid.  Coincidence  of  the  vectors  means  no  /5- 
turning.  Two  numbers  given  at  each  grid  indicate  the  maxi¬ 
mum  /5-turning  between  two  consecutive  levels  such  as 
from  370-m  (level  6)  to  635-m  (level  7)  in  the  sub-tropical 
regions  (10°-30°N)  and  from  1025-m  (level  8)  to  1575-m 
(level  9)  in  the  mid-  and  high-latitudes  (north  of  40°N).  At 
some  locations,  the  maximum  turning  occurs  at  deeper 
levels.  Usually,  strong  vertical  turning  of  the  horizontal 
velocity  indicates  a  strong  baroclinicity. 


70  m,  but  weaker  in  the  P-vector  inverse  method.  The  lower 
level  (1575  m,  and  3250  m)  velocity  fields  are  very  weak. 
The  major  difference  between  the  MOM  and  the  inverse 
solutions  is  the  western  boundary  flow  bifurcation.  The 
MOM  model  shows  the  high  latitude  westward  flow  bifur¬ 
cated  at  40°N  for  the  depth  of  1,575  m,  and  at  60°N  for  the 
depth  of  3,250  m  into  two  western  boundary  currents  (north¬ 
ward  and  southward).  However,  the  P-vector  inverse  model 
does  not  show  this  bifurcation. 


8.  Absolute  Velocities  Obtained  from  the  P-Vector 

Inverse  Method  Using  MOM  T,  S  Output 

Taking  the  statistically  steady  potential  temperature 
and  salinity  fields  as  no-error  data  sets,  we  employed  the  P- 
vector  inverse  method  to  obtain  the  3-D  absolute  geostrophic 
velocity  (ur,  v/).  Figure  11  shows  the  vector  plots  of  the 
horizontal  velocities  at  several  different  depths:  (a)  12.5  m, 
(b)  37.5  m,  (c)  70  m,  (d)  215  m,  (e)  1575  m,  (f)  3250  m.  The 
circulation  patterns  are  very  similar  to  the  MOM  statistically 
steady-state  velocity  fields  in  the  upper  four  levels  (12.5  m, 
37.5  m,  70  m,  and  215  m).  The  circulation  patterns  can  be 
outlined  as  follows.  Westward-moving  equatorial  currents 
are  evident  with  a  width  of  6°  in  latitude.  When  the  equato¬ 
rial  current  approaches  the  western  boundary,  it  turns  direc¬ 
tion  and  become  the  western  boundary  current.  The  current 
velocities  are  similar  at  the  three  levels:  12.5  m,  37.5  m,  and 


9.  Comparison  between  Non-Divergent  Portions  of 
Two  Flow  Fields 

Given  that  the  MOM  model  flow  field  derives  from 
primitive  equation  dynamics,  we  don’t  expect  that  the  ve¬ 
locity  fields  from  the  P-vector  method  (assuming  geostro¬ 
phic  dynamics)  should  match  the  MOM  model  velocity 
fields.  The  only  valid  comparison  should  be  with  the 
nondivergent  portion  of  the  flow  field,  i.e., 


V2  _  MOM  )  djju  mom) 

V  r  MOM  -  -) 

dx 

vV/  =  ^/)  d^fu^ 


dy 


dx 


dy 


(27) 
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Fig.  11.  Absolute  velocity  fields  computed  by  the  P-vector  method  at  different  depths:  (a)  12.5  m,  (b)  37.5  m,  (c)  70  m,  (d)  215  m, 
(e)  1,575  m,  and  (f)  3,250  m. 


We  solved  the  two  Poisson  Equations  in  (27)  by  taking 
(mmom,  vmom),  (ui,  vi)  as  forcing  functions  under  the  bound¬ 
ary  conditions 

V^MOM  |r  =  V/l\r=Q  (28) 

where  T  represents  the  lateral  boundaries,  and  n  denotes  the 
outgoing  normal  direction. 

The  i//mom  field  has  the  following  features  (Fig.  12):  a 
nearly  basin- wide  anticyclonic  gyre  in  the  upper  levels  (12.5 
m,  37.5  m,  70  m,  and  215  m)  and  a  nearly  basin-wide 
cyclonic  gyre  in  the  lower  levels  (1575  m,  3250  m).  These 
gyres  are  evident  by  closed  x/Zmom  contours  with  a  maximum 


value  for  the  anticyclonic  gyre  and  with  a  minimum  value 
for  the  cyclonic  gyre.  Both  basin-wide  gyres  are  asymmet¬ 
ric.  The  center  of  the  gyres  (both  upper  and  lower  levels)  is 
towards  the  west,  with  a  strong  western  boundary  current.  In 
the  upper  levels,  a  weak  and  narrow  cyclonic  gyre  appears 
in  the  high  latitudes  (north  of  60°N)  east  of  55°W  (south  of 
Greenland  and  Iceland.)  This  weak  high  latitude  cyclonic 
gyre  reduces  its  size  with  depth  and  becomes  an  anticyclonic 
gyre  which  stretches  westward  to  the  western  boundary  in 
the  intermediate  level  (1,575  m).  In  the  deep  level  (3,250  m), 
the  cyclonic  gyre  fills  the  whole  basin. 

The  streamfunction  I///  (Fig.  13)  at  the  six  different 
depths  shows  that  the  pattern  of  the  two  fields  are  quite 
similar:  a  nearly  basin- wide  anticyclonic  gyre  in  the  upper 


Evaluation  of  P-Vector  Method  195 


MOM  Psi  Function 


Fig.  12.  Horizontal  distributions  of  ty^MOM  at  different  depths:  (a)  12.5  m,  (b)  37.5  m,  (c)  70  m,  (d)  215  m,  (e)  1,575  m,  and 
(f)  3,250  m. 


levels  (12.5  m,  37.5  m,  70  m,  and  215  m)  and  a  nearly  basin¬ 
wide  cyclonic  gyre  in  the  lower  levels  (1575  m,  3250  m). 
These  gyres  are  featured  by  closed  i ///  contours  with  a 
maximum  value  for  the  anticyclonic  gyre  and  with  a  mini¬ 
mum  value  for  the  cyclonic  gyre.  The  center  of  the  basin¬ 
wide  gyres  (both  upper  and  lower  levels)  is  towards  the  west, 
and  has  a  strong  western  boundary  current.  The  patterns  are 
very  similar  between  i/Zmom  and  if//  at  all  depths.  The  differ¬ 
ence  between  i//mom  and  y/i  are:  (a)  the  nearly  basin-wide 
gyre  is  stronger  in  the  MOM  model  than  in  this  P-vector 
model  except  for  the  intermediate  level  (1,575  m),  where  the 
cyclonic  gyre  obtained  by  the  {3-  vector  model  is  stronger;  (b) 
in  the  upper  levels  (12.5  m,  37.5  m,  70  m,  and  215  m)  the 
high  latitude  (north  of  60°N)  cyclonic  gyre  obtained  by  the 


P-vector  model  stretches  from  the  eastern  boundary  to  the 
western  boundary;  (c)  in  the  intermediate  level  (1,575  m), 
the  high  latitude  anticyclonic  gyre  by  the  P-vector  model  is 
located  in  the  eastern  part  (east  of  40°  W)  rather  than  stretching 
to  the  western  boundary  in  the  MOM  model. 

Relative  difference  (I¥)  between  two  streamfunctions 
(V^mom  and  if/f)  for  each  level  can  be  depicted  by  a  ratio 
between  root-mean- square  (RMS)  difference  and  standard 
deviation  of  MOM  results  (°VM0M), 
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Fig.  13.  Horizontal  distributions  of  y/j  at  different  depths:  (a)  12.5 


Here  Mis  the  total  number  of  the  horizontal  grid  points.  The 
numerator  and  denominator  of  (29)  show  the  mean  differ¬ 
ence  between  two  streamfunctions  and  the  variability  of 
Ymom,  respectively.  The  smaller  the  Iv,  the  smaller  the  dif¬ 
ference  between  the  inverse  method  and  MOM  solutions. 
Near  surface,  Iw  is  small  (-0.15),  increases  with  depth  until 
1000  m  (level  8)  with  the  maximum  value  of  0.8,  and  then 
decrease  with  depth  (Fig.  14).  In  deep  levels  (z  <  -2,500  m), 
I¥  reduces  with  depth  from  0.3  to  0.2.  Near  bottom,  Iw  is 
around  0.2.  The  maximum  value  of  Iv may  caused  by  strong 
baroclinicity  near  that  level  (Fig.  9). 


L'UrpjilLHta 

fl) 

m,  (b)  37.5  m,  (c)  70  m,  (d)  215  m,  (e)  1,575  m,  and  (f)  3,250  m. 


10.  Conclusions 

(1)  Two  necessary  conditions  were  discussed  in  this 
paper  for  the  validity  of  any  inverse  method.  They  are  (a) 
non-coincidence  of  potential  density  and  potential  vorticity 
surfaces;  and  (b)  /5-turning.  The  P- vector  concept  provides 
a  logical  way  of  checking  these  two  conditions.  Existence  of 
the  P-Vector  guarantees  the  satisfaction  of  the  first  neces¬ 
sary  condition.  Existence  of  vertical  turning  in  the  horizon¬ 
tal  P-Vector  components,  Ph  =  ( Px ,  Py),  guarantees  the  sat¬ 
isfaction  of  the  second  necessary  condition. 

(2)  The  P-vector  inverse  method  with  pre-required 
conditions  was  evaluated  using  the  MOM  model.  The  statis¬ 
tically  steady  solutions  of  temperature  and  salinity  from 
MOM  are  used  as  a  no-error  data  set  for  computing  absolute 
geostrophic  velocities  by  the  P-vector  inverse  method. 
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Fig.  14.  Vertical  variation  of  the  relative  difference  Iy. 


Similarity  of  i/Jmom  and  i/// fields  is  found  at  different  depths 
and  confirms  that  the  inverse  method  has  capability  of 
picking  up  the  major  signal  of  the  velocity  field. 

(3)  The  introduction  of  P- vector  makes  the  {3- spiral 
formulation  analytical  and  simple.  However,  the  P-vector  is 
calculated  by  differentiation  of  potential  vorticity  which 
may  suffer  from  various  noises.  The  continuation  of  this 
study  should  include  the  {3- spiral  method  and  compare  among 
P-vector  method,  /5-spiral  method,  and  MOM  results. 
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